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Abstract 



The purpose of this paper is to present a theoretic and numerical study of utilizing squeezing and phase 
shift in coherent feedback control of linear quantum optical systems. A quadrature representation with 
i- G ■ built-in phase shifters is proposed for such systems. Fundamental structural characterizations of linear 

quantum optical systems are derived in terms of the new quadrature representation. These results reveal 
considerable insights of issue of physical realizability of such quantum systems. The problem of coherent 
quantum LQG feedback control studied in [331 150] is re-investigated in depth. Firstly, the optimization 
U" 1 methods in [331 150] are extended to a multi-step optimization algorithm which utilizes ideal squeezers. 

Secondly, a two-stage optimization approach is proposed on the basis of controller parametrization. Nu- 
merical studies show that closed-loop systems designed via the second approach may offer LQG control 
performance even better than that when the closed-loop systems are in the vacuum state. When ideal 
C*") ' squeezers in a close-loop system are replaced by (more realistic) degenerate parametric amplifiers, a suf- 

m 

of performance convergence is also discussed in the LQG control setting. 



ficient condition is derived for the asymptotic stability of the resultant new closed-loop system; the issue 
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1 Introduction 

Being an important branch of quantum physics, quantum optics has been undergoing an accelerating growth 
in its applications to emerging quantum technology as it offers building blocks for constructing quantum 
computing, communication and metrology devices to realize the dream of quantum nano-scale technology 
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[H132]. As sound approximations to the fundamental field models in quantum optics, linear quantum optical 
systems have been developed in terms of quantum stochastic differential equations (QSDEs) [23] (here the 
word 'linear' refers to the linearity of the Heisenberg equations of motion of system operators), based on 
which a vast body of measurement-based feedback control methods have been proposed to achieve such 
objectives as entanglement preservation, state preparation, error correction, etc. Interested reader may 
refer to [3j [12] H6j HH HH 07] for excellent discussions on measurement-based feedback control and their 
applications to a broad range of quantum optical systems. 

An alternative feedback mechanism, coherent feedback control, has been proposed recently where mea- 
surement is not necessarily involved in feedback loops; instead, quantum information may flow directionally 
as a (possibly non-commutative) signal (such as a quantum optical electromagnetic field or an injected laser), 
or directly via a bidirectional physical coupling. The benefits of coherent feedback include (i) preservation 
of quantum correlation of the whole network, and (ii) high speed |244 Fig. 1]. In [45] an all-optical feedback 
mechanism is studied for a quantum optical system comprised of two cavities (a source cavity and a driven 
cavity) ; By designing appropriate interaction Hamiltonian coupling, the authors were able to obtain squeezed 
state inside the source cavity. In [49] Yanagisawa and Kimura derived closed-loop linear quantum system 
models consisting of cavities and beamsplitters. A scheme is proposed in [40] to produce continuous-wave 
fields or pulses of polarization-squeezed light via coherent feedback. This proposal is further investigated in 
|37| . A more general quantum modeling framework is studied in [20], which is Markovian when channel-to- 
channel time delay is ignored. This modeling framework is further studied in [18] where transfer functions 
are obtained for quantum feedback optical networks mediated by beamsplitters. A Hamiltonian formulation 
of such modeling framework is proposed in [19J. Lately the linear quantum systems framework studied 
in [29] and [19] has been extended to include squeezing components [22]. Input-output maps and trans- 
fer functions are defined in this more general framework. In [25] quantum H°° control of linear quantum 
stochastic systems is developed, where the resulting controllers can be classical, fully quantum or quan- 
tum/classical mixed. Based on the theoretical work developed in [25], a quantum optical experiment is 
implemented recently [29]. Because a quantum controller itself is a quantum system, its time evolution 
must obey Schrodinger's equations. To deal with this fundamental issue, the concept of physical realizability 
is proposed in [2S]. Physical realizability is also investigated in [35], [30], [3S], [5T)] . 

As one of the major methods of modern control, linear quadratic Gaussian (LQG) feedback control 
has been extended into the quantum domain. In the framework of measurement-based feedback control, 
quantum LQG feedback control has been investigated in [121 ttU S3. USl UOl SH] , [SI Sec. 6], etc.. Recently, 
this problem has been studied in the setting of coherent quantum feedback networks [33J and [SD]. It turns 
out that coherent LQG control is more challenging than coherent H°° control in that control performance 
and physical realizability are not separable in the LQG coherent feedback control setting. The resulting 
optimization problem is in general non-convex, and analytical solutions are therefore very difficult to find, 
if not impossible. A numerical procedure is proposed in [33] which shows that there exists a fully quantum 
linear controller which offers better closed-loop LQG control performance than classical ones do. A similar 
procedure is proposed in [50] which utilizes direct couplings between plants and controllers. Generally 
speaking, direct coupling is a physical mechanism by which a quantum plant and its quantum controller can 
exchange energy directly, without necessarily through field connections (cf. |26l Sec. II-C and Fig. 4], [5U|, 
Sec. H-B]). Direct coupling can be implemented via nonlinear optical devices like crystals [451 Fig. 1]. 
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An ideal squeezer can be modeled as a static Bogoliubov transform [221 Sec. II-C], which is an idealization 
of a (more realistic) degenerate parametric amplifier (DPA). Ideal squeezers have been used for theoretic 
study in quantum optics [T6l [211 E21 [22] • Unfortunately there is as yet no rigorous theoretical justification 
for their use. Since squeezing devices are prevalent in quantum optics, it is important to address this issue. 

Phase modulators are optical elements that manipulate optical waves at both classical and quantum 
levels; for the latter, they manifest as effects on creation operators of optical modes. Phase shift has been 
used in a variety of quantum applications, e.g., [271 1211 [Hj [36] . 

The purpose of this paper is to study quantum optical networks including squeezing components and 
phase modulators from a control theoretical perspective, aiming at providing systematic control techniques 
for a wide spectrum of applications. The follows three paragraphs outline the major contributions of the 
paper. 

Firstly, we present a quadrature representation of linear quantum optical systems with built-in phase 
shifters (Sec. 12. 5p . This representation contains the usual amplitude-phase quadrature representation as a 
special case. Fundamental algebraic characterizations of such quantum systems are presented in terms of 
this new representation (Theorems [TJ and [2]) . These results reveal more insights of the concept physical 
realizability of open quantum systems first explored in [25]. Theorems Q] and [2]) are the theoretical basis of 
the subsequent numerical investigation. 

Secondly, quantum LQG coherent feedback control is re-studied in the general framework of linear 
quantum optical systems including squeezing components and phase shifters presented in Sec. [2j Firstly, 
we generalize the numerical procedures in [331 [50] by including ideal squeezers and direct couplings in 
closed-loop systems, and show that performance can be improved considerably (Sec. 14. 2p . Secondly, we 
propose a controller parametrization approach and a two-stage optimization technique to find a coherent 
feedback controller, ideal squeezers, direct coupling, and/or phase shifters simultaneously. It is shown that 
appropriate co-design of ideal squeezers, direct coupling, and/or phase shifters can build a closed-loop system 
which offers considerably good control performance (Sees. 14.31 and I4.4p . 

Finally, a sufficient condition is derived for the stability of closed-loop systems when ideal squeezers are 
replaced by DPAs (Theorem[3|). Moreover, a case study is conducted in Sec. !5.2l to demonstrate performance 
convergence in the LQG feedback control setting. This study hopefully will shorten the gap between the 
existing theoretical results and their applications, thereby providing experimentalists with some guidance 
for the choice of parameters of DPAs. 

The rest of the paper is organized as follows. In section [2j the class of linear quantum optical systems 
of interest is presented. Section [3] presents the set-up of closed-loop systems. Section H] studies coherent 
quantum LQG control. Section [5] studies performance of closed- loop quantum systems when ideal squeezers 
are replaced by DPAs. The paper is concluded by Section [6l 

Finally, some words for notation. 

Notation. Let i = y/—l be the imaginary unit. Given a column vector of operators x = [ x\ ■■■ x m ] T 
where m is a positive integer, define = [ x\ • • • x* ] T where the asterisk * indicates Hilbert space ad- 
joint or complex conjugation. Furthermore, define the doubled-up column vector to be x = [ x T (x#) T ] T ■ 
The matrix case can be defined analogously. The symbol diag n (M) is a block diagonal matrix where the 
square matrix M appears n times as a diagonal block. Given two matrices U, V € C rxfc , a doubled-up 
matrix A(U,V) is defined as A(U,V) := [ U V; V# U# ]. Let I n be an identity matrix. Define 
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J n = [ I n ; —I n ] and \I> n = diag(/ n , — /„). (The subscript "n" is always omitted.) Then for a matrix 
X e C 2nx2m , define X b := V m X^ n . Finally, the norm symbol || • || represents the largest singular value for 
a constant matrix. 

2 Linear quantum optical systems with squeezing components 

This section introduces systems of interest in the paper. A quadrature representation with embedded phase 
shifters is introduced, some fundamental relations are presented, an example is used to demonstrate that 
phase modulation is useful in measurement-based feedback control of linear quantum optical systems. 

2.1 Boson fields 

An m-channel free-space light field is described by a vector of annihilation operators 6j n (i) = [bi n ,i{t), • • • , &m,m(£)] r 
defined on a Fock space [IB], |5Uj . These operators satisfy the singular commutation relations 

[binj(t), b* n , k (t')] = 5 ]k 5{t - t'), [bi^t), b in>k {t')} = [b* mJ (t), bt njk (t')) = 0. 

The operator bi n j(t) (j = 1, . . . , m) may be regarded as a quantum stochastic process; in the case where 
the field is in the vacuum state, this process is quantum white noise. The integrated process Bi n j(t) = 
Jq bin t j(r)dT is a quantum Wiener process [HI Sec. 5.3.5], which may be used to define quantum stochastic 
integrals, with associated non-zero Ito product dBi n j{t)dB* nk {t) = 5j k dt. In this paper we assume that 
there is no scattering among channels. 

2.2 Open quantum optical systems 

An open quantum optical system G is a collection of n quantum harmonic oscillators a(t) = [a\(t), . . . , a n (t)] T 
(defined on a Hilbert space) interacting with m boson fields bi n (t). Such system can be parameterized by 
a triple (I m , L, H). In this triple, the vector operator L is defined as L = C^a + C + o*, where C_ and 
C+ E fcrrixn anc j a = . . . , a n ] T with a.j being the annihilation operator of the jth quantum harmonic 
oscillator satisfying the canonical commutation relations [aj,a£] = Sjk- Define a matrix C = A(C_,C+). 
The observable H = iatA(ft_,f) + ) a is the initial internal energy of the oscillators, where f2_ and £7+ & 
C nxn satisfy = Ol and f2 + = JlT. With these parameters, G on the composite system-field space can 
be written as, in Ito form, 

da{t) = -(iVH + ^C l 'C)d(t)dt-C l >dBi n (t) , o(0) = a, 
dB out (t) = Ca{t)dt + dB m (t). (1) 

2.3 Ideal squeezers 

In quantum optical experiments, boson fields bi n (t) can be squeezed before they interact with quantum 
optical systems. Squeezing components can be approximated by Bogoliubov transformations. Let S- and 
S + G C mxm . Define S = A S+). If S satisfies S S = SS b = l2 m , then S is called a static Bogoliubov 
transformation. This paper focuses on a particular class of static Bogoliubov transformations. Let S be a 
real diagonal matrix. Then it can be verified that 5 = exp (A (0, £)) is a static Bogoliubov transformation. 
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In this paper, such type of Bogoliubov transformations are called ideal squeezers, which have been used in 
the study of quantum optics |144 |16j. Clearly, S is an identity matrix if £ is a zero matrix, that is, there is 
no squeezing. 

Example 1 An ideal squeezer can be regarded as an idealization of a phase-shifted degenerate parametric 
amplifier (DPA). A DPA is an open oscillator that is able to produce squeezed output field. A model of a 
DPA is as follows (pH page 220]): 



da(t) 
dB out (t) 



a(t)dt — \fridBi n {t), 



Ka{t)dt + dBi n (t), 



where k and e are assumed to be real numbers for simplicity and satisfy the inequality |e| < k. In terms of 
(I,L,H) language, 0_ = 0, tt + = f , C_ = 0c, and C+ = 0. Let r = lnf±f. Then S = exp(A(0,r)) = 
A (cosh r, sinh r) is an ideal squeezer. 

Replace k and e by £ and j- respectively where h > 0, a sequence of DPAs is obtained parameterized 
by h. It can be shown by simple algebra that the sequence of DPA (parameterized by h) converges to e lw S 
pointwisely as h — > 0. 



2.4 Open quantum optical systems with ideal squeezers 

If the input fields bi n (t) pass through ideal squeezers S before interacting with a collection of open quantum 
harmonic oscillators, then the composite system reads 



da(t) = Aa(t)dt + BSdB in (t) , o(0) = a, 
dB 0Ut (t) = Ca(t)dt + SdB in (t), 



(2) 



in which the coefficient matrices are given by 

1 



A 



c'c -m n A(n^,n + ) , b = -c\ 



(3) 



2.5 Quadrature representation of quantum optical systems with embedded phase shifters 

So far, we have used the annihilation and creation operators a,j and a?j to represent oscillators systems, via the 
doubled-up notation a = [a T a)] T . There is an alternative form, amplitude-phase quadrature representation, 
in which all the operators are observables (self-adjoint operators) and all the corresponding matrices are 
real, not imaginary. 

Phase shift is a technique which can be used to change phase of a light beam. Phase modulators are 
widely used in experimental quantum optics [2]. In this paper we show that systematic design of phase 
shifters can help to improve closed-loop control performance considerably. In this subsection, we introduce 
a new type of quadrature representation with built-in phase shifters. 

Let P n be a 2n x 2n permutation matrix which transform a column vector d = [d±, . . . , d2 n ] T to P n d = 



[g?i, d%, e?5, . . . , d2,d±, . . . d2 n ] T ■ Let be a unitary matrix defined via 



M :-- 



1 

V2 



1 1 








, M k := M 




— i i 




e 







1 

72 



e 
-ie 



-i6 k 
-iOk 
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e> s > 









<b) 



Figure 1: (a) quadrature representation; (b) quadrature representation with embedded phase shifters 

where 6k is real. Introduce 

M a = diag (Mi, ...,M n ),M b = diag (M n+1 , M n+m ) , M c = diag (M n+m+1 , . . . , M n+2m ) . 
Denote matrices 

A a = P re M Pj, A b = P m M h Pl, A c = P m M c Pl. (4) 
It can be readily verified that AtA,- = A,-At = I,j = a, 6, c. Define a coordinate transform 

x := A a a, B in ■= A b B in , y := A c B out . (5) 

Then in quadrature form, the system G in Eq. ((2J) is converted to 

dx(t) = Ax(t)dt + BSdB in (t), x(0) = x, (6) 

dy{t) = Cx{t)dt + A c A\SdB in {t), (7) 

where 

A = A a AA\, B = A a BA\, C = A C CA\, S = A b SAl (8) 

Remark 1 When all 6^ are 0, the above quadrature transforms reduce to the unitary transforms used in 
[501 Eq. (24)] 0. 

The above representation contains the amplitude-phase quadrature representation as a special case. In 
fact, it is the amplitude-phase quadrature representation of the system from bi n to b ou t in Fig. QJb), not that 
in Fig.[Qa). That is, input, output and intra-cavity fields are implicitly assumed to be possibly phase-shifted 
in this new quadrature representation. This new quadrature representation enables us to choose suitable 
quadratures which may yield desired (closed-loop) control performance. 



1 In quantum optics, given annihilation and creation operators a and a* , amplitude operator is denned to be a + a 1 while 
phase quadrature —i(a — a^), cf. [2] Sec. 4.3.1], [421 Eqs. (2.56), (5.5), (5.6)]. Therefore, the quadratures defined by Eq. ([5]) 
are the scaled versions of commonly used quadratures (with the scaling factor l/y/2). However, it can be easily seen that these 
two types of transformations give rise to the same real matrices. 
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Example 2 The following linear quantum plant is studied in |33l Sec. 8]: 



dx(t) 



A 

-A 



+ 




-2y/iti 

2^ 


z(t) = x(t) + (3 u (t), 



dy(t) 



x{t)dt + 

dB in ,i{t) + 
x(t)dt + dB inA (t) 



du(t) 








-IJni 



dBin^it), 



(9) 



where A = 0.1, k\ = k 2 = K3 = 0.01, (3 U is the signal part of u (cf. [25j Eq. (3)]). The quantum LQG 
control problem studied in [33] is to design a controller that minimizes the following performance index: 



1 /"* 

J 00 =limsup- / (z T (t)z(t)) dr. 



(10) 



That is, the aim of control is to steer the quadrature operators as close as possible to the origin in phase-space 
x with a minimum controlling force (5 U for energy consideration. 
Define a quadrature representation 



x = A a a, u = AbB in , B iri) i = AB irij i, B iriy2 = AB iri)2 , y = A c B out} i, 



where 



Aa 



1 


e i9i 


e -iex 


, At 


1 


e i02 


V2 


-ie idl 




"71 


_ ie id 2 


1 


1 1 


, A c = 


1 


e ie 3 


e -i0s 


V2 


—i i 


V2 


_ ie ie 3 


ie -i6s 



-102 
-102 



Then, in quadrature representation the system is 



dx = 

y = 

z = 

where j3 u is the signal part of u, and 



Axdt + Bdu + BidBi ni i + B 2 dBi n<2 , 
C 2 xdt + D 2 idB in} i, 



(11) 



A 



Bl = 2y/Kl 



A 
-A 



sin(0i) 
-005(6*1 



, B = 2^ 

C 2 = 2^ 



- sin(6 ) i) sin(^2) sin cos (#2) 
cos(6'i) sin(6 ) 2) — cos(6 ) i) cos (#2) 



005(6*3) cos(^i) cos (O3) sin (6>i ) 
sin(6 ) 3) cos(6 ) i) sin^) sin(6'i) 



D 



21 



008(6*3) -sin(0 3 ) 
sin(0 3 ) 005(6*3) 



Z = 1,2. 
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To compare with results in [331 Sec. 8.2], only the first element of y is measured. When 9\ = 62 = 63 = 0, 
the system reduces to Eq. (38) in [33]. It has been shown [331 Sec. 8.2] that a measurement-based feedback 
controller yields a closed-loop LQG index 4.8468. In what follows we study several other cases. 

(i) . When 9\ = 62 = and #3 7^ 0, we have continuous measurement by homodyne detection [331 [2]. 
Numerical study shows that 63 = gives rise to an LQG performance Jqo = 4.8468. That is, there is no 
improvement. 

(ii) . When 6± = 0, 2 = -0.5294, and 63 = -0.5498, Joo = 4.0551, which is not only better than 4.8468, 
but also better than 4.1793 — the best coherent LQG control performance obtained in [33]. 

(iii) . When B\ = 0.5200 and 6*2 = 03 = 0, 3Joo = .7544. 

(iv) . When 6>i = 0.04, 6 2 = -0.49, and 63 = -0.1, = 3.7388. 

This example demonstrates that measurement-based feedback plus phase shift can yield LQG control 
performance better than that via coherent feedback studied in [33] and |50| . 

2.6 Fundamental Relations 

In this subsection characterizations of open quantum systems are presented. 

The following result reveals a fundamental structural relation that the quantum systems under study 
have. 

Theorem 1 The real matrices A, B, and C in Eq. satisfy the following equations: 



AJ n + J n A T + BJ m B T = 



(12) 



B = -U n C T (k c m m k\). 



(13) 



Proof. It is easy to show that matrices A and C in Eq. (|3|) satisfy 



y n A + A^ n + C^ m C = 0, 



which leads to 




(14) 



Jn-A + A J n + C Jm,C — 0. 



(15) 



Moreover, by B = -& we have k a Bt^ h = A ( A c CAi) {A c $> m Al). Consequently 



B = -iJ n C T (A c ^ m Al), 



(16) 



which is exactly Eq. (fT3|) . Moreover, by Eq. (fT6|) . C 
above into Eq. (fT5|) gives 



iA c ty m A^ b B T J n . Substitution of B and C derive 



J n A + A T 




J n A + A T J n - J n B A h m m A\ J m A c ^ m A\B T J ? 



n 




= 0. 



(17) 
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Hence Eq. (fT2"j) holds. The proof is completed. □ 

The following result characterizes solutions to the equations (|12p - (|13p . This result connects equations 
(fl~2j) - (fT3]) with open quantum systems of the form (pQ). 



Theorem 2 Given matrices A a ,Afc, A c defined in Eq. there is a unique matric C of the form C = 

A(C_,C+) and infinitely many matrices A of the form A = —^C^C — if n A (fi_, 0+) such that matrices 
A,B,C defined via Eq. ||}) satisfy Eqs. (T ^) -(T ^) . 

Proof. We first show that the matrix A has the specified form. By Eqs. (I12D and (]17p we have J n A + 
A T J n + C T J m C = 0. As a result, the matrix A is in the form of A = ^J n C'J m C + J n R, where R 



is 



an arbitrary real symmetric matrix. Consequently, in terms of A = AaAA a and C = aJca a we have 

a . However it is straightforward to show that a is in the form of A (fi_ , 

Finally, since R is arbitrary, there are infinitely many A of the form A = — ^C^C — i^f n A Next we 

show that the matrix C has the specified form. Given a real matrix i?, some simple algebraic manipulation 
gives that -iJ n BA b ^ m = \X* X] for some complex matrix X. Substituting it into Eq. (|13p yields 
C = A C [X X*] T . Consequently, 



c = k\CK 



X T 



A a = A(C_,C+) 



for some matrices C_ and C+. The proof is completed. □ 

Remark 2 Eqs (|12p and (|13p are forms of physical realizability investigated in [25] which reveals fun- 
damental relations that a quantum system must obey. The derivation of physical realizability in [25] is 
much mathematically involved. Interested reader may check the proofs of Theorems 2.1 and 3.4 in [25]. 
In contrast, the derivation presented here is much simpler as we make use of connections between the 
annihilation-creation and quadrature representations. Moreover, Theorem [1] extends the results in [25] to 
linear quantum optical systems with built-in phase shifters. Therefore, the results presented here are more 
general, cf. [251 Remark 3.6]. These results are used in the numerical study carried out in Sec. [H 

Remark 3 The indication of Theorem [2] is two-fold. Firstly, for a given triple (A a ,A&,A c ), the matrix 
C is uniquely determined by Eqs. (|12p and (|13p . while A has free parameters. That is, an open quantum 
system of the form (pQ) solving Eqs. (|12p and (|13p has unique input-output coupling structure while its initial 
internal Hamiltonian is allowed to be arbitrary. Secondly, the free choice of triple (A OJ A&, A c ) provides us 
with freedom of choosing quantum systems that solve Eqs. (|12p and (|13p . This advantage is utilized in the 
numerical study in Sec. 14.41 



3 Closed-loop systems 

This section presents the set-up of the closed- loop systems. 

Consider the closed-loop system as shown in Fig. [2l where P is a quantum plant and K a coherent 
feedback controller to be designed. S u , S y , S v ki, S v k2 are ideal squeezers to be designed. 6j n> i and 6^2 are 
quantum noises, bi n 1 is in vacuum state, while bi n ^ may have finite energy, y is plant output (output field 
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b,„, 



H SyK1 h 



byk10 



Figure 2: Coherent feedback control 



channel), u is control input (input field channel), u a is part of the output of the controller K. z is reference 
output (namely, a performance variable, which may not be a physical variable, so is designated by a dotted 
line), b Vkl and b Vk2 are quantum vacuum noise inputs to the controller K. Hi nt in Fig. [2] designates the 
direct coupling between P and K, cf. [50l Eq. (7)]. 

In terms of quadrature representation introduced in Sec. I2.51 the quantum plant P is described by 



dx(t) 
dz(t) 
dy{t) 



and the controller by 

dec*) 

du (t) 



= Ax(t)dt + Bdfi(t) + BtdBin^t) + B 2 dB iri)2 (t), x(0) = x, 

= C lX {t)dt + Du(t) + D l2 dB in ^(t), 

= C 2 x{t)dt + D 2 idBin,i(t) + D 22 dB in)2 (t), 

A K £(t)dt + B K dy {t) + B K1 dB VKW (t) + B K2 dB VK20 {t), £(0) = & 
C K i{t)dt + b K dB VK10 {t). 



(18) 



(19) 



Note that in the quadrature representation introduced in Sec. 
as can been seen from Eq. ([7]) and the case study in Sec. H7 
Define matrices 



Dk in general is not an identity matrix, 



A BS u Ck + B\2 

BkS v C 2 + B 2 i Ak 



B d 



B 2 
BkS v D 22 



(20) 



G, 



cl 



H, 



d 



B\ BS u DkS v ki 
B K S y D 2 i Bk\S v ki Bk 2 S v k2 

o bs u b K s vK i o 



C\ DS u Ck 



where B\ 2 and B 2 \ are for direct coupling, satisfying B 2 \ = J n Bf 2 J n with B\ 2 being an arbitrary real-valued 
matrix [50] . The closed- loop system with a built-in direct coupling in the quadrature representation is given 
by 



dx{t) 



dz{t) 



A, 



Cr 



x(t) 

m 

x(t) 



dt + B d dB inj2 (t) + G c 



dt + D 12 dB in o(t) + H c 



dBi n ^{t) 

dB VK1 (t) 
dB VK2 (t) _ 

dB in ^(t) 
dB VK1 (t) 
. dB VK2 (t) 



(21) 
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4 LQG synthesis 

In the previous section we have added ideal squeezers and phase shifters to closed-loop plant-controller 
systems. It is natural to investigate if they are helpful for the purpose of control. As mentioned in the 
Introduction section, coherent LQG quantum feedback control is as yet an outstanding problem, in this 
section we would like to re-study this problem in order to illustrate the usefulness of the more general 
framework presented in the previous sections. Example [2] has shown that measurement-based feedback plus 
phase shifters can outperform the best coherent feedback control performances derived in [33] and [50]. We 
show in this section that coherent feedback can achieve much better performance if ideal squeezers, phase 
modulators and direct coupling are appropriately designed in coherent feedback control. 



4.1 Set-up 



In this subsection, we formulate the coherent quantum LQG feedback control problem. This problem has 
been investigated in [33] and [50]. Similar problems have been investigated in [12], [11], [33], [15], [TO], [38] , 
\A7\ Sec. 6], etc., in the framework of measurement-based feedback control. As mentioned in Example [21 the 
aim of control is to minimize the variances of amplitude and phase quadratures with a minimum controlling 
force for energy consideration. 

The following assumptions are standard (cf. 



Al. D\2 in the plant (|18p is zero. 

A2. The quantum noise inputs 6m,i(i) and bi n ^{t) in the plant (|18p are in vacuum state. 

A3. The plant is initially in Gaussian state. 

A4. The controller is initially in vacuum state. 

Moreover, for simplicity, we assume 

A5. The performance variable is z{t) = x(t) + /3 Uo (t), where f3 Uo (t) is the signal part of u {t). 

Under these assumptions, suppose that A c \ in Eq. (|20p is Hurwitz, then a positive definite matrix Vlqg 
is the (unique) solution to the following Lyapunov equation 



A c iV L q G + V LQG A T cl + 



B c l G c l 



B c l Gel 



0. 



As is shown in [33], the closed-loop LQG performance index in (|10p can be computed via 



Joo = Tr[C d V LQG C ( 



T 



(22) 



(23) 



As we wish to minimize the variances of amplitude and phase quadratures, to conform to footnote 1 in 
Sec. 12. 51 we use the following quadrature transformation 



(24) 



1 1 




a 


—i i 




a* 
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In what follows we first study an extreme case. Denote 

Cl C2 1- (25) 

C2 C 3 

Theorem 3 Assume that in steady state both the plant and the controller are in the vacuum state, then 
Joo = 2 + (ci +c 3 ). 

Proof. Firstly, it is easy to show that 

(z T (s)z(s)) = (x T (s)x(s)) + (x T (s)C K as)) + (e(s)C^x(s)) + (f{s)C%C K ti{s)) . (26) 

Secondly, note that 

x T (s)x(s) = 4a* (s)a(s) + 2, (27) 

and 

(e(sf C K C K ^s)) = (c iq l(s)+c 2 q^s) P ^s) +p^ s )q^s) + c 3 pl(s)) , (28) 

where 

gf = a^ + l + 2a* Y OK + (a* Y ) 2 , p£ = - (a^ - 1 - 2a* K a K + {a* K ) 2 ) , 

qm+Pm = - 2i (a 2 K - (a* K ) 2 ). 

Substituting Eqs. ([27]) and ([28]) into Eq. {26} yields 

<5 T (s)z(s)> = -2ic 2 (a 2 K (s) - ( a * K (s)) 2 ) + (c, - c 3 ) (a 2 K (s) + (a* K (s)) 2 ) (29) 

+ (x T (s)C K as)) + (e(s) T C-i:x( S )) + 2 + ( Cl + c 3 ). 

As both the plant and the controller are in the vacuum state in steady state, 

<a( S )*o(*)> = 0, (a 2 K (s)) = ((a* K (s)) 2 ) = 0, (x T (s)C K t(s)) = (s)C],x(s)) = 0. (30) 

Substituting Eq. ([30]) into Eq. (J29]) leads to (z T (s)z(s)) = 2 + (ci + c 3 ). Consequently, = 2 + (ci + c 3 ). 
The proof is completed. □ 

Remark 4 It is worth pointing out that derived above is a result of the Heisenberg's uncertainty 
principle. Clearly, is independent of plant parameters. 

Remark 5 The plant studied in Example [2] can model an optical cavity interacting with three input fields. 
If the term /3 U is ignored in Eq. ([9]) by setting Ck = 0, then it is easy to see that the performance index 
defined in Eq. (|10|) means minimization of the number of photons in the cavity. Theorem shows that the 
minimal value is 2. However, the plant in this example is only marginally stable, so it has to be controlled. 
Therefore Ck cannot be zero and hence = 2 + {c\ + c 3 ) > 2. 

As direct coupling can add new energy channel between the plant and its controller, squeezers can modify 
the input-output structure of the plant, and phase modulators can adjust relative phases of signals in the 
closed loop, in this sequel we seek to study the following problem. 

Problem 1 Is it possible to design a coherent feedback controller, direct coupling, ideal squeezers and phase 
modulators simultaneously such that the resulting closed-loop system performance better than 2 + (ci + c 3 )? 

The measurement-based feedback studied in Example [2] fails to provide an affirmative answer to the 
above problem. Moreover, the coherent controllers proposed in [33] and [50] fail too. 



C k Ck = 



12 



4.2 An Extension of algorithms in 



In this subsection we re-investigate the example studied in |33| 150] . It turns out that an appropriate 
choice of ideal squeezers and direct couplings can improve the performance of the fully quantum controller 
considerably. This case study embodies the essential ingredients of quantum LQG controller design, thus 
algorithms for more general cases can be developed without additional conceptual difficulty 
In the following we present a numerical optimization algorithm. 

Let 2 by 2 real matrices Zj 1( . . . , Zj i2 1; Zj 14 1 be a set of decision variables matrices. Denote 

V = [J; Zj i; ■ ■ ■ ; Zj i2 x ; Z^ x ; ■ ■ ■ ; Z^ 14 x ] T . Define a symmetric matrix Z = VV T . For ease of presentation, 
denote 



Zx 2 ,i Z X3 ,i Z X4 ^i 
SvKl 

o S vK2 



5 D W 







5„ 








// 







, D w = 




. <?2 . 





B2 BS u S v ki 

SvKl 

5'„_ft'2 







S„ 
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. 




. An 









Construct the following linear matrix inequalities 



AZ Xlhl +BS u Z Xs , 1 +(AZ a . atl +BS u Z XB ,i) T Z^i+A B w 

Z Xl ,i+A T Z X7A A + Z x cl + (Z XrA A + Z x c'i) T Z X7 sB w + Z x D n 

(B'J T (Z X7tl B' w +Z x Dj T -I 



+ 



a rank constraint 



B12 + Bf 2 
Z xl12 ,i+Z X7> iBi2 




Z Xl2t i+Z X7j iBi2 

+ 




Z xil ,i + z xiil 



<o, 



Z X6 ,1 

I 



C\Z x&) i + Di2Z X6j i C\ 



I \C\Zx 6 ,\ + Di2Z X5t i 
Z X7 ,i &i 

Q 

Tr(Q) < 7, 
rank(Z) < 2, 



>0, 



(31) 

(32) 
(33) 
(34) 
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and a set of additional constraints 



Z > 0, 



Z X gA + Z x 



0. 



z 



0,0 



l,x 6 



Z Vg A ~ 

Zvu,l 
Zvi2,l 

Zvi,l - 

Z X g,l ~ 

Zv 5 ,l ~ 
Zx 12 ,l 



1 = 0, 

- Z X6 ,i = 

- Z X7 ,i = 

~ Zxi^XQ " 
Zv4 ,VQ 

Zv$ ,V-J 



Z Vi ,i — Z X7 aB^Su — 0, 



0. 
0, 

= 0, 
= 0, 
= 0, 
= 0, 



z 



"6,1 



z, 
z , 



Zv\,X2 ~ 

Z X2 ,iJn v 



+ Z Vg> l 



0. 



0, 
0. 
0. 

= 0, 

Zv2,X3 — 0, 
Zv3, X4 = 0, 

Zx 3 ,iJn Vk2 = 
Zx4,iJn v = 0, 



Zy 7 ,l *-*Xg,XQ 

Zv%,i Z Xl X g - 
Z Xl0 ,i + Z\ )Xl0 
Z Vl3 ,i 

Zvi4,l 
Zv2,l ' 
Zv3,l - 



(35) 



Z XA ,\C2 — Z X7y \A — 0, 



z x 



xn,l 



Z X g,\B2\- 



Two equality constraints for physical realizability: 



z V8 ,i + zj s l + z Vlli i - zj iul + Z Vl0t i 



Z Vl0 ,i + Z Vl2 ,i + Z Vri x + Z Vlij i 



and 



Zx-2,1 



Z Vr aJn v 



0. 



(36) 
(37) 



The quantum LQG control algorithm is given below. 



Initialization. Set B\2 = 0, S u = I, S v = I, S, 



VKl 



I, and 5, 



t'A'2 



Step 1. For fixed B12, S u , S y , S VK1 , and S VK2 , given 7 > 0, employ a semidefinite programming (SDP) 
approach to solve the feasibility problem with constraints Eqs. (|3ip - (|37|) in which decision variables are 
Z and Q. 

Step 2. Pertaining to Step 1, use some optimization procedure to find direct coupling parameters B\2- 
Step 3. Pertaining to Step 2, use some optimization procedure to find squeezers S u , S VK1 and S VK2 - 
Then go to Step 1. 

Step 1 can be solved using an algorithm similar to that proposed in [33] which is based on LMIRank 
([55]). SeDuMi ([38]) > and Yalmip ([28]), while in Steps 2 and 3 only a few decision variables are involved, 
hence they can be easily solved via many general-purpose optimization algorithms. 

If the above optimization problem is solvable, controller parameters of K can be obtained via Ak = 
Z Xl ,i, Bki = Z X2t i, Bk2 = Z X3t i, Bk = Z X4t i, Ck = Z X5i i. 

In what follows we present the solution of the preceding numerical algorithm. After initialization, solving 
step 1 yields a controller K with parameters 



.4 



K 



0.0251 
0.0665 



BK2 



10 



-10 



-0.3787 
-0.2121 

0.6008 
0.1938 



Br 



1.0273 
0.8235 



-0.1964 
-0.0492 



) Bki 



0.1125 
0.1504 



-0.5992 
-0.1284 



-0.3049 
-0.2273 



Ck 



0.1284 
0.1506 



-0.5993 
-0.1126 



, D K = I. 
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which yields = 4.1787. Implement step 2 produces a direct coupling with 



B 



12 



1.2 -9 
0.72 0.36 



x 10 



-3 



B 



21 



-0.36 
0.72 



-9 
■1.2 



x 10" 



The it is found that Joo = 3.9995. Implementing step 3 generates squeezers 



S U Sy /, Sy Kl 



1.5876 
0.6299 



VK2 



1.8076 
0.5532 



It can be verified that in this case = 3.8312. This is, by adding appropriate squeezers LQG performance 
is further improved. With direct coupling and squeezers obtained above, run step 1 again to yield a controller 
K with parameters 



A K = 
Bri 
C k = 



0.0108 -0.4819 

0.0353 -0.2116 
0.2996 -0.7662 
0.1427 -0.1876 

0.1876 -0.7662 

0.1427 -0.2996 



Br - 

U.i 

B K 2 = lO" 14 
Dk = L 



1.3697 -0.4995 
0.7232 -0.1559 

0.9902 0.6591 
0.4085 -0.0584 



and the resulting LQG performance is = 3.7464. 



4.3 Controller parametrization and a numerical optimization approach - without phase 
shifters 

The algorithm in the preceding subsection is an extension of those in [331 I50| . however such algorithms 
suffer from severe limitation: they are nonlinear and non-convex optimization procedures involving many 
decision variables. Consequently, it is very challenging to use such algorithms to optimize over all controller 
parameters, direct couplings, ideal squeezers, and phase shifts simultaneously. 

From Sec. 12. 21 we see that controllers can be parameterized by C_, C+, f2_, 0+, from Sec. 12 .31 we see that 
each ideal squeezer can be parametrized by a real number; and direct coupling and phase shifters involve 
just a couple of parameters. Moreover, with such parametrization physical realizability is naturally satisfied. 
This enlightens us to seek for controller design methods via optimization over these parameters. Inspired 
by the study in [7] , in this and the next sections we investigate a new optimization approach on the basis of 
controller parametrization and a two-stage optimization algorithm. 

The idea of controller parametrization is simple. The coherent feedback controller to be constructed is in 
the form of Eq. (|19p . that is, a quantum harmonic oscillator interacting with three input fields. According 
to Sec. 12.21 there are 13 parameters to be optimized. Direct coupling requires 4 parameters, while each 
of the four squeezers needs one parameter respectively. Therefore, we will optimize over 21 parameters 
simultaneously. Based on this parametrization, a two-stage optimization is constructed which is outlined as 
follows. 

In Stage one, a genetic algorithm optimizes over the cost function (|23p subject to the constraint that 
the closed- loop matrix A c i is Hurwitz. 
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Note that physical realizability is satisfied naturally in terms of parametrization. The genetic algorithm 
performs a global search over the parameter space to find a minimal solution |17j . 

In Stage two, a local search is performed around the minimal solution derived from Stage one. Around 
this solution, the closed-loop matrix A c i is always Hurwitz, hence the optimization problem boils down to 
minimizing the cost function (|23|) subject to the constraint (|22|) . Here, the cost function and constraints are 
all continuously differentiable. Consequently, this optimization problem can be solved by means of many 
standard optimization algorithms, cf. [31]. 



Remark 6 In contract to the numerical optimization procedure presented in UcK [hU\/ and the previous sub- 
section, the above optimization problem has a parameter space of much smaller dimension. Moreover, genetic 
algorithms usually provide a reasonably good initial solution 

The preceding numerical algorithm turns out to be very effective. It provides the following system 
parameters 



+ 



-0.0136 + 0.0857i 
-0.0473 - 0.3509i 
2.7099 + 19.4445i 
n_ = 0.9768, n + = -2.4874 - 0.3771?', 

Therefore, the matrices of controller K are 



0.0136 + 0.0857i 
0.0286 - 0.225H 
0.0763 - 0.3437i 



A K 



-193.0685 
1.5106 



3.4642 
192.3144 



0.0272 0.0000 
0.1715 0.0000 





0.1715 



, Bk - 

Br2 = 



-0.0272 



-2.6336 
19.1008 



-19.7882 
-2.7862 



0.0759 0.1258 
-0.5760 0.0187 

b K = i. 



The matrices for direct coupling are 
B\2 = 

Finally ideal squeezers are 

Si i 



151.1269 0.0621 
1.3904 123.8024 



B 



21 



230.3001 
0.0043 



-123.8024 
1.3904 



0.0972 




0.0621 
-151.1269 





10.2920 



S 



VK-2 



0.5163 
1.9367 



1.1253 
0.8887 



The resulting LQG performance is 2.0004, which is much better than that obtained in Sec. 14.21 

Finally, by Eq. (|25p . 2 + (c\ + C3) = 2.03013 > 2.0004. That is, this closed-loop system offers performance 
even better than the vacuum case given by Theorem 
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4.4 Controller parametrization and a numerical optimization approach - with phase 
shifters 



In this subsection we study whether the LQG performance can be further improved if phase modulators are 
added. According to the study in the previous section, the controller, direct coupling and ideals squeezers 
obtained there already yields a performance which is better than that in the vacuum case, thus adding phase 
modulators will not improve performance considerably. 

As shown in Example [21 the parametrization of the quantum plant requires three phase values, say 6\, 62 
and #3. As far as the controller is concerned, we denote by #4, #5, and 6q the phase shifters of the intra-cavity 
mode of the controller, the traveling fields y Q and u Q respectively. That is, there are 6 more optimization 
variables than the optimization problem studied in Sec. 14.31 

Using a similar two-stage optimization procedure as that in Sec. 14.31 we find 0\ = 0.2711, 62 = —1.0493 x 
1(T 4 , 6» 3 = -0.30233, = 3.8780, 0+ = 2.7264 + 0.3709i and 



C_ = W 2 x 



0.0010 + 0.0032i 
-0.0178 + 4.5222i 
0.0567 + 0.0117i 



0.0474 + 0.3290i 
1.7123 - 6.7882i 
0.5665 + 0.2253i 



The matrices for direct coupling are 

Bu = 10 5 x 
Finally, the ideal squeezers are 



1.3774 0.002003 
0.005402 0.4648 



B 2 i = 10 5 x 



10 3 x 



9.2358 

0.0000001083 



-0.4648 0.002003 
0.005402 -1.3774 



0.4808 
2.0798 



S 



VK2 



1.7476 
0.5722 



, Sy 



1.5537 
0.6436 



The resulting LQG performance is 2.0000008422, which is slightly better than 2.0004. 



Remark 7 The big difference between parameters of controller, direct coupling and ideal squeezers obtained 
in the above two subsections reveals the non-convex nature of the underlying optimization problem. This in 
turn reveals the necessity of simultaneous optimization. 



5 Performance convergence 

The preceding section has shown that ideal squeezers are very useful in controller design for linear quantum 
optical systems. In the practice of quantum optics, ideal squeezers are approximations of more realistic 
degenerate parametric amplifiers (DPAs). Therefore, it is desirable and practically important to check how 
closed-loop quantum LQG performance changes when ideals squeezers are replaced by DPAs. In this section, 
we focus on closed-loop stability and closed-loop LQG performance convergence. 
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5.1 Asymptotic stability 



The concept of stability is discussed in |5CH Sec. III]. In particular, for the linear quantum systems discussed in 
this paper, asymptotic stability is equivalent to that the "A"-matrix is Hurwitz stable. In this subsection we 
study asymptotic stability of closed-loop systems when ideal squeezers are replaced by degenerate parametric 
amplifiers (DPAs). 

For simplicity, we assume that D22 = in Eq. (|18|) . As in Example Q3 S u is replaced by 



dx u (t) 
du(t) 







x u (t)dt - \/-^du (t) := A u x u (t)dt + B u du (t) 



-jj-x u (t)dt + du (t) ) := C u x u (t)dt + D u du (t). 



(38) 



Assume other ideal squeezers S y , S v ki, S v k2 are replaced in a similar way. Then the closed- loop system 
becomes 



dx 




d£ 








dx v k2 




d>Xyo 









A 

BxD y C2 + B21 



ByC2 



BD u Ck + B12 BD u DKC v k\ 



A K 






BRlCykl 

A v ki 







BK2C v k2 BkC v 





A v k2 





+ 








B u Cr 


BuDxCvkl 





Bi 


B 2 


BD u DxD vkl 







B K D y D 2 i 





BuiDyki 


BR2D v k2 




dBi n ^i 








B v ki 







dBin t 2 











B v k2 




dB vk i 


ByD2l 













dB vk2 








BuDxDvki 








The following matrices are used in the sequel. 



Ai 



A 

BkS v C2 + B21 



BS u Ck + B 



12 



A 



A 



A 2 



ByC2 












B u Ck 



BC U 







-4, 



BtfCy 



BC U 






Ay 









, A 4 = 


A u 


, X = 



VhA 2 A 2 A 3 A 4 1 



Notice that T is independent of h, while ^3 contains 1/h. 



x 

Xykl 
X v k2 



X 



!/<> 



(39) 



Theorem 4 Assume the closed-loop system \21\) is asymptotically stable. Then the closed-loop system \39\ 
is asymptotically stable for all h > if the matrix T + T T is negative definite. 
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Proof. The asymptotic stability of the closed-loop system (|39p is the same as that of the following system 



dx 




A 


BD u Ck + B12 





BC U 




X 






BkD v C 2 + B 2l 


A K 


BrC v 









d>Xy 




B y C 2 





Ay 









dx u 







B u Cr 





A u 







(40) 



Si 


£ 2 


BD u DxD vkl 







dBi ny i 


BkD v D 2 i 





BxiDyki 


BK2D v k2 




dBinZ 


ByD 2l 













dB vk x 








BuDKD vk i 







dB vk2 



+ 



Therefore, in what follows we focus on system (|4U|) . Let us look at its "A" -matrix. Some algebraic manip- 
ulation yields 



7 

/ -B K C y Ay l 

J 





I 







/ B K C v A~ l 








I 










J 

BC U A U 1 


7 



BxD y C 2 + B 2 i 

ByC 2 





BD u Ck + B\ 2 
Ak 


B u Ck 





BrC v 

Ay 





7>C« 





.4 



BS u Ck + Bx 2 (BS u Ck + B^BxCyA 1 



BkS v C 2 + f?2i Ar: 

ByCa 

B u Ck 

Ax AxA-sA^ 1 
A 2 A A + A 2 AzAl l 



A K B K C y A- x 



A, 



BuCRBxCyAy 1 



ABCuA- 1 
(BkS v C 2 + B 2 i)BiC u A~ 
By^BCuA^ 1 
A*, 



Moreover, 

Then by Fact 5.12.3 in @] 
ReA 



/ 




A x AxA 3 A 4 1 




' I 


-i 





Vh 




A 2 A4 + A 2 A 3 A^ 




Vh 




A A 



+ T. 



-i 



Ax 7^1^4 

VhA 2 A± + A 2 A Z A~, X 



— 2 ^max 











( 





2^4 



+ -A max (T + T T ) 



(41) 



Note the right-hand side of (|4ip is independent of h. Consequently, if the matrix T + T T is negative definite, 
then the closed- loop system ([39]) is asymptotically stable for all h > 0, so is the closed- loop system (j2~T|) . □ 
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5.2 Performance convergence 

Pertaining to the study in Sec. 14.31 if we replace the squeezers S u , S v ki, S v k2, and S y by DPAs respectively 
then the closed- loop system (j39[) is obtained. According to Theorem simple calculations show that this 
closed-loop system is asymptotically stable for all h > 0. Fig.[3]shows how the closed-loop LQG performance 
changes as a function of h. It can be seen that, as h goes to zero, the performance converges to 2.02588, 
which is slightly bigger than 2.0004, but still better than 2.03013, the vacuum case given by Theorem [3J 

2.09 1 1 1 1 1 1 

2.08 ■ _ — " - 

2.07 
2.06 

8 

2.05 
2.04 
2.03 



Figure 3: Closed-loop quantum LQG Performance with replacement of ideal squeezers by degenerate para- 
metric amplifiers 

6 Conclusion 

Recent years have seen a considerable amount of work on quantum networks theory, e.g., system intercon- 
nection @9l UHl UHl EQl [22], H°° control [251 ED], LQG control [331 ED], synthesis theory [SHIM], dissipation 
theory and direct coupling [26l [50] , laboratory demonstrations [29], among others. This paper fits into such 
general picture by presenting how to use phase shift and squeezing in control design for linear quantum 
optical feedback networks. In addition, the controller parametrization and two-stage optimization proce- 
dure may be useful in the optimal design of quantum optical networks, their advantage for engineering 
non-classical correlation has been illustrated in the study of coherent LQG quantum feedback control. 
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